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چکیده 


استخراج بیش از حد از آب‌های زیرزمینی می‌تواند منجر به کمبود آب شیرین» نفوذ آب‌شور به آبخوان‌های ساحلی و درنتیجه شوری بیش از حد 
آن‌ها شود. مدل‌های شبیه‌سازی آب زیرزمینی ابزاری مفید جهت شناخت رفتار سیستم آبخوان‌های ساحلی هستند. سواحل جنوبی دریای خزر یکی از 
مناطقی است که با معضل پیشروی آب‌شور به آبخوان‌های ساحلی روبرو است. به‌منظور شبیه‌سازی جریان آب زیرزمینی و پیشروی آب شور در آبخوان 
ساحلی دریای خزر در منطقه ساری- نکا مدل‌های عددی 10۳۳۲01۷ ۷1۲3۳315 و ٩۳۸۷۸‏ در محیط نرم‌افزار 071510.0 به کار 
گرفته شدند. پس از شبیه‌سازی برای ارزیابی مدل, واسنجی در مدل کمی برای تراز آب زیرزمینی و در مدل کیفی برای غلظت شوری به مدت ۴ سال از 
مهر ۱۳۸۹ تا شهریور ۱۳۹۳ به‌صورت ماهانه انحام گرفت. مدل واسنجی شده با استفاده از اطلاعات در دسترس برای یک سال اینده از مهمر ۱۳۹۳ تا 
شهریور ۱۳۹۴ برای تراز آب زیرزمینی و غلظت شوری صحت‌سنجی شد. نتایج پهنه‌بندی مکانی مدل حاکی از کیفیت مطلوب آب زیرزمینی در نواحی 
مرکزی و جنوبی آبخوان ساری-نکا بود. پس از صحت‌سنجی مدل و با فرض ثابت ماندن شرایط هیدروژئولوژیکی آبخوان» نتایج پیش‌بینی مدل برای 
نواحی با 1195 بیشتر از ۲۰۰۰ میلی گرم بر لیتر برای ۶ سال آینده بیانگر هجوم آب‌شور در مناطق شمالی آبخوان بود. 


واژه‌های کلیدی: آبخوان ساحلی. آبخوان ساری-نکاء تهاجم آب‌شور غلظت شوری 


مقدمه 

پیشروی آب‌شور در آبخوان‌های ساحلی با توجه به حجم بالای 
برداشت از منابع آب زیرزمینی در این نواحی و با توجه به قابلیت‌ها و 
پتانسیل‌های بهره‌برداری به‌عنوان یکی از معضلات مدیریت منابع آب 
تبدیل شده است. در نواحی ساحلی آبخوان آب شیرین با آب‌شور 
ارتباط هیدرولیکی داشته و با توجه به میزان تراز آب زیرزمینی جهت 
جریان به سمت دریا یا اقیانوس برقرار است. تغییر در میزان تراز آب 
سبب معکوس شدن جریان آب زیرزمینی می‌شود و اين موضوع سبب 
می‌گردد تا جبهه‌های آب‌شور وارد آبخوان شود. این تهاجم و انتقال 


۱ و ۲- به‌ترتیب دانش آموخته کارشناسی ارشد و استادیار مهندسی آب و سازه‌های 

هیدرولیکی» دانشکده عمران, دانشگاه صنعتی نوشیروانی بابل بابل ایران 

(:- نویسنده مسئول: 1۲ ۵ 1 تححقط :تم 
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آب‌شور از طرف دریا به ساحل به‌عنوان تداخل آب‌شور و شیرین یا 
نفوذ آب‌شور معرفی می‌شود. تهاجم جبهه‌های آب‌شور می‌تواند 
تأثیرات منفی زیست‌محیطی بر سیستم بهره‌برداری ساحلی بگذارد و 
بهره‌برداری از آن‌ها را به‌عنوان منبع تأمين آب با کیفیت مناسب 
به‌منظور تأمین آب شرب و انجام فعالیت‌های کشاورزی محدود کند. 
با توجه به افزایش جمعیت و افزایش نیاز از یک‌طرف و تنش‌های 
اقلیمی و محیطی از طرف دیگر مدیریت منابع آب شیرین آبخوان‌ها 
بسیار حائز اهمیت است؛ بنابراین با توجه به مشکلات تهاجم 
جبهه‌های آب‌شور در نوار ساحلی» توجه به مدیریت آبخوان‌های 
ساحلی جهت کاستن از حرکت گوه آب‌شور به آبخوان‌ها و بالازدگی 
مخروطی " آب دریانزدیک چاه‌های پمپاژبهعنوان یکی از ولویت‌های 
اصلی مدیریت منابع آب مطرح است. 


صمتمتاتاصاً 52110۳۷26 -4 
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۰ نشریه آب و خاک. جلد ۰۳۴ شماره ۲. خرداد - تیر ۱۳۹۹ 


درجه نفوذ آب‌شور در مکان‌ها و شرایط هیدروژئولوژیکی مختلفه 
متفاوت است. در بسیاری از موارد» ناحیه آلوده‌شده توسط نفوذ آب‌شور 
به بخش‌های کوچکی از آبخوان محدود می‌شود و هیچ تآثیری بر 
روی چاه‌های پمپاژ که جهت تأمین آب مورد استفاده قرار می‌گیرد 
نداشته و یا تآثیر بسیار کمی دارد. در سایر موارده آلودگی به‌صورت 
منطقه‌ای بوده و تأثیر قابل‌توجهی بر منابع آب زیرزمینی دارد. 
گسترش آب‌شور به یک آبخوان به چندین عامل بستگی دارد. این 
عوامل شامل مجموع آب زیرزمینی تخلیه‌شده نسبت به مجموع آب 
رین بشید ایام فصو له تس ها تهاسان انا لها 
زهکشی) از منبع! (منابع) آب‌شور ساختار زمین‌شناسی و توزبع 
خصوصیات هیدرولیکی آبخوان و حضور واحدهای محدودکننده‌ای که 
ممکن است مانع حرکت عمودی آب‌شور به سمت آبخوان یا درون 
آبخوان شوند. می‌باشد. علاوه بر اين» زمان لازم برای جبهه آب‌شور 
جهت عبور از آبخوان و رسیدن به یک چاه پمپاژ می‌تواند بسیار 
طولانی باشد. این زمان بسته به محل و عرض جانبی ناحیه انتقالی 
دارد (۱۶). هجوم آب‌شور استفاده از آب‌های زیرزمینی را جهت 
مصارف شرب کشاورزی پا صنعت محدود می‌کند. ازاین‌رو ضرورت 
پیش‌بینی مکان و حرکت آلودگی در برابر خطرات احتمالی وجود دارد 
(۱۰). طی سه دهه اخیرء» پژوهش‌های متعددی در ارتباط با جنبه‌های 
مختلف حرکت املاح در آبخوان‌ها صورت گرفته است. همچنین 
مدل‌سازی پیشروی آب‌شور به سیستم آب زیرزمینی توجه زیادی را به 
خود اختصاص داده است و چندین مدل ریاضی و عددی در اين ارتباط 
توسعه یافته است. این مدل‌ها به‌منظور پیش‌بینی مرز مشترک يا ناحیه 
انتقالی بین آب زیرزمینی (آب شیرین) و آب دریا در مناطق ساحلی به 
کاز گرفته می‌شوند: (۱۹) توعد این مل‌ها عمدتا به دلینل مشافل 
مختلف مربوط به آب‌های زیرزمینی, ازجمله ارزیابی ذخایر آب 
زیرزمینی» پیش‌بینی نفوذ آب‌شور و حرکت روبه‌جلو یا رو به بالای مرز 
مشترک در پاسخ به شیوه‌های بهره‌برداری از آب زیرزمینی بوده است. 

مسئله پیشروی آب‌شور در سراسر جهان ازجمله در مادراس واقع 
در هند (۲۳). آبخوان دلتای نیل در مصر (۲۴) حوزه آبریز یون لین در 
جنوب غربی تایوان (۲۷) مناطق جنوب غربی بنگلادش (۲۰)» هماندو 
در فلوریدا 4)٩(‏ حوزه آبریز رودخانه جائه در شاندوگ چین (۴) و ساوانا 
در گرجستان (۱۴) گزارش شده است. در بیشتر این مطالعات عامل 
اصلی پیشروی آب‌شور به آبخوان‌های ساحلی» بهره‌برداری بیش‌ازحد 
و عدم مدیریت صحبح منابع آب زیرزمینی بوده است. لازمه مدیریت 
آبران‌های ساحلی زاشخاب استزاتی‌هآی مناسب یت کشرل و 
کاهش نفوذ آب‌شور می‌باشد. این استراتژی‌ها تنها زمانی می‌توانند 
تکامل یابند که فرآیندهای فیزیکی موجود در آبخوان‌های ساحلی 
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به‌طور قابل‌اعتمادی شبیه‌سازی شده و مورد مطالعه قرار گیرند. 
مدل‌های ریاضی به کار گرفته شده برای یک منطقه خاص پس از 
واسنجی و صحت‌سنجی می‌تواند جهت شبیه‌سازی پاسخ سیستم 
آبخوان با در نظر گرفتن استراتژی‌های مختلف پمپاژ و تغذیه به کار 
گرفته شوند. در سالیان اخیر محققان بسیاری به بررسی نحوه پیشروی 
در رها و ای بای شاه از سل ها ناف 
شبیه‌سازی پرداخته‌اند که از آن جمله می‌توان به مطالعات زیر اشاره 
داشت. 

وفایی و همکاران (۲۶) با به‌کارگیری مدل 5۳۸۷۷۸۲ به 
بررسی چگونگی پیشروی آب‌شور در اثر برداشت بی‌رویه آب زیرزمینی 
در سواحل پرداختند. اعمال هفت نرخ افزایشی مختلف برای تراز آب 
نشان داد که اين مسئله کاهش گرادیان هیدرولیکی آب زیرزمینی و 
افزایش نفوذ شوری به آبخوان ساحلی را به دنبال دارد. همچنین انجام 
پمپاژ با دبی دو و سه برابر و احداث مانع زیرزمینی به ترتیب منجر به 
افزايش و کاهش پیشروی شوری شده است. کتابچی و همکاران (۱۵) 
در پژوهش خود با به‌کارگیری مدل عددی ٩۳115۸‏ به بررسی 
تأثیرات تغییرات اقلیم در میزان پیشروی آب شور در آبخوان‌های 
ساحلی شیب‌دار پرداختند. نتایج نشان می‌دهد که تأثیر تومان افزایش 
سطح آب دریا و تغییرات خط ساحلی موجب تشدید پیشروی آب شور 
می‌گردد که در مناطق ساحلی با شیب کم» پیشروی آب شور بیشتر 
است. رضاپور و همکاران (۲۲) در تحقیق خود چگونگی پیشروی آب 
شور را تحت شرایط کاهش سطح ایستابی آبخوان مورد بررسی قرار 
دادند. نتایج بررسی آن‌ها نشان داد که دو عامل اختلاف چگالی بین 
آب شور و شیرین و شیب هیدرولیکی معکوس شده سبب پیشروی آب 
شور به داخل آبخوان می‌گردد. فتحی زاد و همکاران (۷) به‌وسیله مدل 
۳ به بررسی موقعیت مرز مشترک آب شور و شیرین در زیر 
حوزه تالار پرداختند. نتایج پژوهش آن‌ها نشان داد که در مناطق غربی 
زیر حوزه تالار آب شور تا فاصله ۴/۵ کیلومتری از خط ساحلی به 
داخل آبخوان نفوذ کرده است. همچنین در نواحی مرکزی میزان 
پیشروی آب شور در حدود ۲/۵ کیلومتر از خط ساحلی می‌باشد. 
کاردان مقدم و همکاران (۱۳) در پژوهشی به بررسی اثرات 
زیست‌محیطی هجوم جبهه‌های آب‌شور به آبخوان سرایان واقع در 
استان خراسان جنوبی پرداختند. نتایج این پژوهش نشان داد که 
مدیریت برداشت از آب زیرزمینی با استفاده از سناریو مدیریتی کاهش 
برداشت آب از منابع زیرزمینی به‌منظور بهبود کیفیت آب آبخوان‌های 
کویری و ممانعت از تهاجم آب‌شور کویری به آن‌ها با توجه به 
خشکسالی‌های اخیر ضروری است. کوبانر و همکاران (۶) در مطالعه 
خود به توسعه یک مدل برای کنترل و جلوگیری از نفوذ آب‌شور در 
آبخوان ساحلی دشت ۹111116-0016 با به‌کارگیری مدل 
7 در طول دریای مدیترانه در ترکیه پرداختند. نتایج 
شبیه‌سازی نشان داد که آبخوان دشت یادشده به افزايش بهره‌برداری 


از آب زیرزمینی حساس است. چیتراکار و همکاران (۵) به بررسی 
الگوی نفوذ آب دربا در دشست ساحلی طهعناه ۸1۳ واقع در عمان 
پرداختند. مدل ساخته‌شده برای شبیه‌سازی دو سناریوی کاهش پمپاژ 
و تغذیه مصنوعی با پمپاژ ثابت برای کاهش نفوذ آب‌شور به کار گرفته 
شد. پیش‌بینی‌های مدل نشان داد که بهینه‌سازی نرخ پمپاژ و تغذیه 
مصنوعی موجب کاهش نفوذ آب دریا در دشت ساحلی طفسناه ۸13 
می‌شود. 

با مروری بر مطالعات صورت گرفته می‌توان به آهمیت بررسی 
مسئله تهاجم شوری در آبخوان‌های ساحلی به‌منظور مدیریت و 
برنامه‌ریزی صحیح جهت کاهش اثرات زیان‌بار آن پی برد؛ بنابراین با 
توجه به اينکه در آبخوان‌های سواحل شمالی کشور آب زیرزمینی 
به‌عنوان یک منبع اصلی تأمین آب شرب و کشاورزی به شمار می‌رود 
و به علت مجاورت این آبخوان‌ها با دربای خزر امکان تهاجم شوری 
به این مناطق ساحلی؛ نابودی منابع آب زیرزمینی و اثرات مخرب 
زیست‌محیطی و اقتصادی وجود خواهد داشت؛ لذا در پژوهش حاضر با 
به کارگیری مدل عددی ۷5) اقدام به ارزیابی کمی و کیفی منابع 
آب زیرزمینی آبخوان ساحلی ساری- نکا واقع در استان مازندران» 
شبیه‌سازی غلظت کل جامدات محلول " و کلراید تعیین مرز مشسترک 
آب شور و شیرین و پیش‌بینی نحوه پیشروی آب‌شور طی زمان در 
آبخوان ذکر شده پرداخته می‌شود. 


مواد و روش‌ها 
منطقه مورد مطالعه 

منطقه مورد مطالعه این تحقیق در محدوده ساری- نکا واقع 
می‌باشد که بین طول‌های جغرافیایی "۵۲۹۲۴ تا ۵۴۹۴۴۳ شرقی و 
۵۶ تا ۲۶۲۵۲ شمالی در شمال کشور ايران قرار دارد. مساحت 
این محدوده ۶۹۳۸/۵ کیلومترمربع می‌باشد که ٩۷۰‏ کیلومترمربع آن 
را دشت و مابقی اراضی در ارتفاعات است. سه رودخانه اصلی تجن. 
نکارود و سیاهرود آبخوان ساحلی را تغذیه می‌کند. روند عمومی 
منحنی‌های تراز آب زیرزمینی در این دشت شرقی- غربی بوده و 
فواصل منحنی‌های هم‌تراز سطح آب زیرزمینی در اراضی مخروط 
افکنه‌ای و جنوبی دشت کم بوده و شیب هیدرولیکی آب زیرزمینی 
بالا می‌باشد. با حرکت به سمت شمال در اراضی میانی دشت و پایاب 
دشتی فواصل منحنی‌ها افزایش یافته و شیب هیدرولیکی آب 
زیرزمینی کاهش می‌یابد. مقادیر توزیع سالانه بارندگی» تبخیر و دما 
طی دوره آماری ۴۵ ساله (۱۳۴۵-۹۰) به ترتیب برابر ۶۶۴ میلی‌متره 
۰ میلی‌متر و ۱۶/۵۱ سانتی گراد می‌باشد (۱۸). موقعیت 
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شبیه‌سازی پیشروی آب‌شور در آبخوان‌های ساحلی ۲۷۱ 


محدوده مطالعاتی در شکل ۱ نشان داده شده است. 
روش پژوهش 

در اين پژوهش ابتدا با به‌کارگیری مدل ۷0۲۲10۷۷ به 
ارزیابی کمی منابع آب زیرزمینی پرداخته می‌شود و نتایج حاصل 
زآن جهت شبیه‌سازی انتقال شوری با در نظر گرفتن فرآیند انتقال, 
پراکندگی. انتشار وارد مدل ۱۷1131218 می‌شود و درنهایت 
خروجی‌های حاصل از مدل 1۷۲13121۷15 جهت شبیه‌سازی انتقال 
شوری با در نظر گرفتن تغییرات چگالی وارد مدل 5۳۸1/۸7 
می‌شود. 


معادلات حاکم 

در این پژوهش از مدل عددی ٩۳۸۷۷۸۲‏ که توانایی 
شبیه‌سازی سه‌بعدی جریان آب زیرزمینی با در نظر گرفتن چگالی 
جریان در محبط متخلخل را دارد استفاده شده است. مفهوم اساسی 
این مدل عددی, ترکیب دو نرم‌افزار شبیه‌ساز جربان آب زیرزمینی 
۷ (۱۱) و شبیه‌ساز انتقال آلاینده 13145 (۲۹) به 
یک برنامه واحد جهت حل معادلات جریان آب زیرزمینی وابسته به 
چگالی و انتقال آلاینده می‌باشد. 10111017 یک مدل تفاضل 
محدود سه‌بعدی برای شبیه‌سازی جربان آب زیرزمینی در شرایط 
ماندگار و غیر ماندگار می‌باشد. معادله حاکم برای جریان آب زیرزمینی 
تحت تآثیر چگالی که توسط ۷1010۷۷ در مدل 5۳۸۷۸۲ 
حل می‌گردد به شرح زیر بیان می‌شود (۸): 
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که در آن 0 م۵ و 7 محورهای مختصات عمود برهم در 
جهت‌های اصلی نفوذپذیری» ,یک روک و مک هدایت هیدرولیکی در 
سه جهت اصلی (۷5)» 0 چگالی سیال ( 4162/0 :۵ چگالی آب 
شیرین ( 2/۳ :هد آب شیرین (0 2 ارتفاع بالای مرکز 
سلول مدل (50 ٩:‏ ضریب ذخیره (1/0)» و تخلخل موثر (بدون 
بعد» ) غلظت املاح (تلعع» م چگالی آب ورودی از منبع تغذیه 
یا خروجی از طریق منبع تخلیه (8/507 ,9 نرخ حجمی جریان بر 
واحد حجم با توجه به تغذیه و تخلیه آب در واحد زمان (/1) و ) زمان 
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5 ساختاری مشسابه مدل ۷0۳۳۲/0۷۷ دارد و 
به‌منظور شبیه‌سازی تغییرات غلظت آلاینده‌های محلول با در نظر 
گرفتن فرآیند انتقال. پراکندگی, انتشار و واکنش‌های شیمیایی به کار 
می‌رود. بعد از توسعه و واسنجی مدل جریان, اطلاعات موردنیاز برای 
شبیه‌سازی انتقال آلودگی (مانند حجم جریان انتقال‌یافته بین سلول‌ها) 
به این مدل داده می‌شود (۲۱). کد ٩۳۸۷۲۷۸۲‏ با به کارگیری مدل 
عددی 1۷۲۲31715 به حل معادلات حاکم بر انتقال املاح که 
ب‌صورت زیر بیان می‌شود» می‌پردازد (۲۸): 
2-۳۸۳۵۲۷66 
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که در آن (1 ضریب پراکندگی هی درودینامیکی (9/) بت 
غلظت املاح ورودی از منبع تغذیه یا خروجی از طریق منبع تخلیه 
(تصلع) و (لا ,... ,1) بط نرخ تولید یا تجزیه در واکنش > از 
واکنش متفاوت (10*5 /ع)» هستند. 

در شبیه‌سازی جریان وابسته به چگالی و انتقال املاح فرض 
می‌شود که چگالی سیال تنها تابعی از غلظت املاح است و از اثرات 
تغییرات فشار و دما بر چگالی سیال صرف‌نظر می‌گردد (۱۷). معادله 
خطی جهت تبدیل غلظت املاح به چگالی سیال در کد 5۳۸۷۷۸۲ 
به‌صورت رابطه (۲) بیان می‌شود: 
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0 
مقدار ج در رابطه فوق بستگی به واحد انتخاب‌شده توسط 


کاربر طی شبیه‌سازی دارد. برای نمونه اگر واحد متر و کیلوگرم برای 
0 
شبیه‌سازی استفاده شود مقدار ِ برابر ۳ ,+ است که 


ب‌صورت تقریبی معادل تغییرات چگالی سیال تقسیم بر غلظت املاح 
برای آب شیرین و آب‌شور می‌باشد. 


مدل مفهومی و مدل‌سازی کمی آبخوان 

مدل مفهومی جربان آب زیرزمینی یک چارچوب توصیفی و 
یکپارچه است که می‌توان با استفاده از آن داده‌ها و اطلاعات مربوط 
به هیدرولوژی آب‌های زیرزمینی را تفسیر نمود. تهیه مدل مفهومی 
مستلزم تلفیق دقیق داده‌هاء اطلاعات و گزارش‌های مربوط به آبخوان 
و جریان آب زیرزمینی در محدوده مطالعاتی است (۲۱). 

پارامترهای لازم برای مدل مفهومی شامل تراز آب چاه‌های 
مشاهده‌ای» موقعیت و دبی چاه‌های بهره‌برداری» دبی رودخانه‌هاء 
توپوگرافی سطح زمین» سنگ کف محدوده مدل‌سازی و هندسه 
آبخوان است. تخمین ضرایب هیدرودینامیکی آبخوان (هدایت 
هیدرولیکی و آبدهی ویژه برای مدل کمی و پراکندگی طولی برای 
مدل کیفی) و تعیین میزان جریان‌های ورودی و خروجی آب زیرزمینی 
و تعبین هرگونه منابع تغذیه و تخلیه آبخوان از دیگر پارامترهای 
موردنیاز جهت تهیه مدل مفههومی است. 


پس از تهیه مدل مفهومی, برای انجام مدل‌سازی عددی شبکه 
تفاضل محدود با سلول‌هایی به ابعاد ۵۰۰ * ۵۰۰ متر مورد استفاده 
قرار گرفت. در گام بعدی هندسه آبخوان که شامل مرزهای مدل و 
توپوگرافی سنگ کف 9 سطح آبخوان است با استفاده از روش‌های 
درون‌یابی در مدل اعمال گردید. لایه‌های تراز ارتفاعی سطح زمین از 
داده‌های 1171 استخراج و با استفاده از درون‌یابی» کلیه نقاط به‌طور 
مستقیم به شبکه مدل تعمیم داده شد. لایه ارتفاعی سنگ‌بستر 
آبخوان نیز به کمک تفاضل داده‌های تراز سطح زمین و نقشه تغییرات 
ضخامت آبخوان ساخته شده و از طریق درون‌یابی وارد مدل گردید. 
مدل‌سازی آبخوان در حالت ماندگار برای مهر ۱۳۸۹ و در حالت غیر 
ماندگار برای یک دوره ۴ ساله شامل ۴۸ دوره تتش از مهر ۱۳۸۹ تا 
شهریور ۱۳۹۳ در نظر گرفته شد. به دلیل کمتر بودن نوسان تراز آب 
زیرزمینی و توازن بین تغذیه و تخلیه ابخوان در ماه اول سال آبی» 
مهرماه برای شروع مدل‌سازی انتخاب گردید. در طی این دوره مدل 
واسنحی شد و صحت‌سنجی نیز برای یک دوره یک‌سلله از مهر 
۳ تا شهریور ۱۳۹۴ انجام گرفت. 

با ترسیم نقشسه هم‌پتانسیل آب زیرزمینی در نرم‌اف زار 
5 با استفاده از آمار تراز آب زیرزمینی ۱۴ حلقه چاه 
مشاهده‌ای در مهرماه ۱۳۸۹ و اعمال آن به سلول‌های مدل از طریق 
درون‌یابی» شرایط اولیه مدل تعیین شد. مرز ورودی و خروجی آبخوان 
شدت جریان وابسته به بار هیدرولیکی 01153 در نظر گرفنه شد. 
بخش‌های جنوبی آبخوان به‌عنوان مرز ورودی و بخش شمالی آبخوان 
که محل تخلیه جریان‌های سطحی رودخانه‌های واقع در محدوده 
مطالعاتی به دریا می‌باشند مرز خروجی جریان در نظر گرفته شد. 
هدایت هیدرولیکی در محدوده آبخوان با توجه به نقشه هم‌قابلیت 
آبخوان و نقشه هم‌ضخامت لایه آبدار به مدل اعمال گردید. با توجه 
به ناهمگنی آبخوان و متفاوت بودن بافت زمین‌شناسی منطقه در نقاط 
مختلف» محدوده مطالعاتی موردنظر به چندین ناحبه جهت اعمال 
هدایت هیدرولیکی و آبدهی ویژه تقسیم‌بندی گردید. 

در شرایط غیر ماند گار» مقادیر آبدهی ویژه براساس نتایج ر 
آزمایش پمپاژ و باتوجه به مطالعات جباری و همکاران (۱۲) برآورد 
شده است. متوسط ضریب هدایت هیدرولیکی در آبخوان ۱۰ متر بر 
روز برآورد گردید. مقدار حداکتره حداقل و متوسط آبدهی ویژه به 
ترتیب ۱۰۸ و ۳/۴ درصد در نظر گرفته شد. مقادیر این دو پارامتر در 
مرحله واسنجی تصحیح خواهند شد. میزان نفوذ از بارندگی در سطح 
ناحیه دشتی ساری- نکا نیز ۶۵/۸۲ میلی‌متر محاسبه شده است که 
۱ درصد از میزان بارندگی سطح دشت را شامل می‌شود (۱۸). 
همچنین با توجه به بافت خاک دشت» وضعیت کرت بندی و شیوه 
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آبیاری ۰ درصد از نرخ برداشت چاه‌های کشاورزی کاهش داده شد و 
تحت عنوان تغذیه از طریق آب برگشتی کشاورزی به مدل اعمال شد 
(۱۸). با توجه به متغیر بودن موّلفه‌های تغذیه و تخلیه آبخوان در 
برای هر دوره تنش متفاوت خواهد بود. 

در دشت ساری- نکا به دلیل غرقابی بودن ابیاری در فصول 
کشاورزی (کشت غالب منطقه برنج است) تبخیر از آب زیرزمینی 
صورت نمی‌گیرد؛ بنابراین در فصول کشاورزی از تبخیر و تعرق 
صرف‌نظر می‌گردد و در سایر فصول مقدار تبخیر از آب زیرزمینی برابر 
۸ میلیون متر مکعب در نظر گرفته شد (۱۸). خصوصیات 
رودخانه‌ها جهت ورود به مدل عددی با توجه به نقشه‌های حد حریم و 
بستر رودخانه‌ها 9 مقادیر آبدهی در گزارشات بیلان منابع آب منطقه 
استخراج و در مدل اعمال شد. در شکل ۲ موقعیت چاه‌های 
بهرهبرداری و مشاهده‌ای کمیء سلول‌های فعال و غیرفسال, مرزها و 
سلول‌های مربوط به رودخانه آورده شده است. 


مدل‌سازی کیفی آبخوان 

پس از ساخت» واسنجی و صحت‌سنجی مدل کمی جریان» مدل 
کیفی آب زیرزمینی در نرمافزار 14731045 تهیه گردید. در مدل 
کیفی که به بررسی پدیده انتقال املاح در آب زیرزمینی می‌پردازد 
مانند مدل جریان زیرزمینی نیاز به تعریف شرایط مرزی است. شرایط 
مرزی در نظر گرفته در این مدل شامل مرز با غلظت اولیه معین که 
در تمام محدوده مدل اعمال شده و مربوط به زمان شروع مدل‌سازی 
است. می‌باشد. مرز با غلظت معین برابر با ۳۵۰۰۰ میلی‌گرم بر لیتر در 
قسمت شمالی آبخوان برای پارامتر 115" و ۲۸۲/۲ میلی اکی‌والان 
بر لیتر (در حدود ۱۰۰۰۰ میلی‌گرم بر لیتر) برای پارامتر کلراید در 
مجاورت دریای خزر است. غلظت اولیه مدل از طریق درون‌یابی 
مقادیر غلظت 11«5 و کلراید قرائت‌شده در چاه‌های مشاهده‌ای کیفی 
(۱۸ حلقه چاه کیفی) در آغاز شیه‌سازی (مهر *۸) به مدل اعمال شده 
است. مقادیر پخشیدگی طولی اولیه. نسبت‌های پراکندگی جانبی افقی 
به پراکندگی طولی و پراکندگی جانین قاقم بسه پراکندگی طولی بد 
ترتیب برابر ۰ متر و ۰/۱ 9 ۰/۱۰۱ در نظر گرفته شده است (۲ و زا" 
مقدار تخلخل در آغاز شبیه‌سازی در تمامی محدوده مدل برابر ۰/۳۵ و 
/ مه ۳۹ جا» / ی 
ضریب انتشار ملکولی برابر ‏ ۱۰ متر مربع بر ثانیه منظور گردید. اين 
مقادیر طی فرآیند واسنحی تصحیح خواهند شد. همچنین اطلاعاتی از 
قبیل محدوده اجرای مدل, شبکه‌بندی و توپوگرافی سطح زمین و 
سنگ‌بستر همانند مدل کمی است. تقسیم‌بندی زمانی کیفی آبخوان با 
گام‌های زمانی ۶ ماهه و با توجه به داده‌های ماه‌های آبان 9 
اردیبهشت در ۸ دوره تنش از مهر ۸٩‏ تا شهریور ٩۳‏ صورت پذیرفت. 
صحت سنجی مدل کیفی به مدت یک سال شامل دو دوره تنش ۶ 
ماهه از مهر ٩۳‏ تا شهریور ٩۳‏ اجرا شد. 
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نتایج واسنجی و صحت‌سنجی مدل کمی 

اصولا یک مدل آب زیرزمینی با عدم قطعیت‌های مختلفی روبه‌رو 
است. یکی از بارزترین پارامترهای ایجادکننده عدم قطعیت در 
شبیه‌سازی» پارامترهای هیدرودینامیکی آبخوان (هدایت هیدرولیکی و 
آبدهی ویژه) می‌باشند که باید به نحوی مطلوب ناحیه‌بندی و مقادیر 
حالت ماندگار برای گام زمانی مهر ۱۳۸۹ با استفاده از داده‌های تراز 
آب مربوط به اين ماه و سپس در حالت غیر ماندگار در دوره زمانی 
مهر ۱۳۸۹ تا شهریور ۱۳۹۲ اجرا و واسنجی شد. این واسنجی با 
استفاده از داده‌های تراز آب مربوط به این دوره آماری به‌منظور به 
حداقل رساندن اختلاف تراز آب مشاهداتی و محاسباتی در محل 
چاه‌های مشاهده‌ای به روش سعی و خطا به ترتیب با تغییر محدود در 
پارامترهای هدایت هیدرولیکی و آبدهی ویژه انجام شد. لازم به ذکر 
است که نتایج مدل واسنجی شده در حالت ماندگار به‌عنوان شرایط 


اشکال ۲ و ۴ پهنه‌بندی مقادیر هدایت هیدرولیکی در واسنجی ماندگار 
و آبدهی ویژه در ابتدا و انتهای واسنجی غیر ماندگار ارائه شده است. 
با توجه به اشکال ۳ و ۴ مشخص می‌گردد مقادیر این پارامترها در 
نواحی شمالی کمتر و در نواحی جنوبی بیشتر است که این امر 
نشان‌دهنده نفوذپذیری کمتر و ریزدانه‌تر بودن اندازه ذرات خاک در 
بخش‌های شمالی آبخوان است. 

مقایسه مقادیر مشاهداتی در مقابل مقادیر محاسباتی تراز آب 
زیرزمینی در دوره واسنجی در شکل ۵ ارائه شده است. نتایج 
به‌دست‌آمده نشان‌دهنده تطابق قابل‌قبول بین مقادیر محاسباتی و 
مشاهداتی است. 

همچنین در شکل ۶ تراز سطح آب مشاهداتی و محاسباتی در 
طول زمان واسنجی غیرماندگار در محل دو چاه مشاهده‌ای به عنوان 
نمونه ارائه شده است. 

مقایسه مقادیر مشاهداتی در مقابل مقادیر محاسباتی تراز آب 
زیرزمینی در دوره صحت‌سنجی در شکل ۷ ارائه شده است. نتایج 
به‌دست‌آمده نشان‌دهنده برازش قابل‌قبول بین مقادیر محاسباتی و 
مشاهداتی است. میزان خطای آماری مدل در مرحله واسنجی و 
صحت سنجی در جدول ۱ آورده شده است. میزان خطاهای آماری 


برآورد شده در شبیه‌سازی مدل کمی حاکی از صحت مدل ساخته‌شده 
می باقن 
ارزیابی مدل کیفی 

فرآیند واسنجی در مدل کیفی همانند مدل کمی به روش سعی و 
خطا و با تغییر پارامترهای پراکندگی افقی» نسبت‌های پراکندگی 
جانبی افقی به پراکندگی طولی و پراکندگی جانبی قائم به پراکندگی 
طولی صورت پذیرفت. در طی مرحله واسنجی مقادیر غلظت 115 و 
گاید قرامت‌شنه در جا‌های نش اهته‌ان مشب با مقادی خاظه 
شبیه‌سازی‌شده توسط مدل طی دوره آماری مهر ۱۳۸۹ تا شهریور 
۳ مقایسه گردید تا میزان خطای آماری مدل به سطح قابل قبولی 
از دقت برسد. جهت کنترل روند فرآیند واسنجی مدل کیفی مقایسه 


شبیه‌سازی پیشروی آب‌شور در آبخوان‌های ساحلی ۲۷۵ 


بین مقادیر غلظت 1195 و کلراید قرائت‌شده در چاه‌های مشاهده‌ای 
قر یل مقادن غاقای ۱5 و کل این تیش فشه رسط میل ور 
شکل ۸ ارائه شده است. همان‌طور که مشاهده می‌شود غلظت آب 
زیرزمینی با ضریب تبیین کمتر از ۰/۸۵ همبستگی کمتری نسبت به 
مقادیر تراز آب زیرزمینی دارند. دلیل این امر تأثیر کم تغییرات 
پارامترهای موّثر بر پخش به‌عنوان تنها عامل ورودی مدل, بر کاهش 
خطای آماری مدل در فرآیند واسنجی است. 

افتکال کر ۱ نب فایسه یم فا ها اسان و 
محاسباتی را در محل دو چاه‌های مشاهده‌ای طی واسنجی مدل کیفی 
نشان می‌دهد. 
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شکل ۷- مقایسه تراز آب مشاهداتی و محاسباتی در صحت‌سنجی مدل جریان 
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شکل ۱۱- مقایسه غلظت مشاهداتی و محاسباتی برای (الف) ۲05 (ب) نآن) بعد از صحت‌سنجی مدل 
۵06 ۲۵۱۱02610 20 ان رظ) ر1۳5 ره 0مامانهاه 0صه ۲60 وهای 0 ممفانهصجص0۵ -11 مناهز۲ 


صحت‌سنجی مدل کیفی جهت تأیید مدل در بازه زمانی مهر 
است. همبستگی بین غلظت‌های 11 و کلراید شبیه‌سازی‌شده در 
محل چاه‌های مشاهداتی و غلظت اندازه‌گیری شده صحرایی در شکل 
۱ ارائه شده است. با توجه به مقادیر ضریب تبیین مقادیر مشاهداتی 
و محاسباتی همبستگی قابل قبولی دارند. 

پس از اجرای نهایی مدل کیفیء برای ارزیابی نتایج مدل از 
شاخص‌های آماری میانگین مجذور خطاء میانگین مطلق خطاو 
میانگین خطا استفاده گردید که به ترتیب برای پارامتر 115 در دوره 
وی ها غترمافکاز یزیا ۳۶ ۳ ۱ میلی کوخ جر رود 
دوره صحت‌سنجی برابر با ٩‏ ۶ ۱۸- میلی گرم بر لیتر محاسبه 
گردید. همچنین مقادیر میانگین مجذور خطاء میانگین مطلق خطا و 


میانگین خطا برهسب میلی‌گرم بر یتر بهترتیب برای پارامتر کلراید 
صحت‌سنحی مدل برابر با ۶۵ ۴۹ ۱۳ محاسبه شد. نتایج مدل کمی 
و کیفی نشان‌دهنده دقت قابل‌قبول شبیه‌سازی است؛ بنبراین از صدل 
می‌توان جچهت پیش‌بینی‌های ی استفاده نمود. 

نتایج حاصل از پهنه‌بندی 115 در شکل ۱۲ ارائه داده شده 
است. نتایج نشان می‌دهد این پارامتر در نواحی نزدیک به ساحل در 
شمال آبخوان دارای بیشترین غلظت و کیفیت نامناسب جهت 
زیرزمینی به سمت جنوب محدوده مطالعاتی کاسته می‌شود. مناطق 


5 کمتر از ۱۰۰۰ میلی‌گرم بر لیتر بخش اعظمی از آبخوان را 


پوشش داده است که اين مناطق حاوی آب زیرزمینی با اندکی شوری 
می‌باشد که برای کشاورزی مناسب هستند. مقادیر ۲15" به سمت 
شمال آبخوان روند افزایشی دارد و در بخش‌های میانی آبخوان که 
عمده چاه‌های بهره‌برداری در این مناطق قرار دارند مقادیر آن از 
۰ میلی گرم بر یتر تا ۲۰۰۰ میلی‌گرم متفیر است. همچنین 
بررسی نقشه مربوط به غلظت کلراید نشان می‌دهد که غلظت این 
متغیر همانند غلظت 115 در بخش‌های شمالی آبخوان و در نواحی 
ساحلی به‌طور چشمگیری افزایش داشته است به‌طوریکه غلظت 
کلراید این مناطق از آب زیرزمینی نسبتاً شور تا آب دریا متغیر است. 

در ادامه به‌منظور مدل‌سازی پیشروی شوری به‌جای در نظر 
گرفتن تمامی محدوده مطالعاتی» مناطتی که غلظت 758 آن‌ها 
بیشتر از مقدار ۲۰۰۰ میلی گرم بر لیتر است در نظر گرفته شده است. 
این انتخاب با توجه به توصیه نشریه ۸۵0 مبنی بر درنظرگیری 
حداکثر 1125 مجاز جهت مصارف کشاورزی انجام شده است. 
شبیه‌سازی با مدل 91۸۱۷۸7۲ و پیش‌بینی وضعیت 
آیخوان 

در اين قسمت به‌منظور شبیه‌سازی انتقال آلودگی تحت تأثیر 


شبیه‌سازی پیشروی آب‌شور در آبخوان‌های ساحلی ۲۷۹ 


تغییرات چگالی از مدل عددی ۳۸۷۷۸۲ در محیط نرم‌اف زار 
0 استفاده شده است. گسسته‌سازی منطقه موردنظر شامل 
شبکه‌هایی با ابتاه ۲۵۰ ۲۵۰ مر در راستای طول وخرضی تما 
بر ۸۰۸ سلول فعال و در ۵ لایه با توجه به خصوصیات منطقه و 
ضخامت آبخوان در نظر گرفته شده است. است. ضخامت آبخوان به 
دلیل وضوح در نمایش ۲۵۰ متر در نظر گرفته شده است. هندسه مدل 
و شرایط مرزی مدل در شکل ۱۳ ارائه شده است. مرز جنوبی آبخوان 
که جریان آب زیرزمینی برای تخلیه به سمت دریا وارد آن می‌گردد و 
همچنین مرز شمالی آبخوان که منطبق بر دریای خزر است (خروجی 
آبخوان) به‌عنوان مرز با جریان وابسته به بار هیدرولیکی در نظر گرفته 
شده‌اند. همچنین غلظت ۲۵۰۰۰ میلی‌گرم بر لیتر برای پارامتر 115 
و غلظت ۲۸۲/۲ میلی اکی والان بر لیتر (حدود ۱۰۰۰۰ میلی‌گرم بر 
لیتر) برای پارامتر کلراید به دلیل مجاورت با دربای خزر به مرز شمالی 
مدل اعمال شد. میزان برداشت از چاه‌های واقع در محدوده در نظر 
گرفته‌شده نیز طبق نتایج آمار استخراج‌شده از شرکت آب منطقه‌ای 
مازندران به مدل اعمال گردید. 
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۵5 6 ()) موتهع۲ 3 (ظ) ونهع۲ 1 (۸) ماه وه ممتاععز ععمتن معط صز بان ۵۲ ممناصطتاعنل ممتاحاصم‌جصم -16 عساهز۲] 


502 رم 

عدد 500 1حو 
502 وحم 

عد 500 1حوح 
500 عم 


شکل ۱۷- تغییر موقعیت خط هم غلظت ۵۰ درصد ب[) در مقطع عرضی 2-9 پس از (الف) یک سال. (ب) سه سال. (ج) شش سال 
6 (0) و۲6۵15 3 (ظ) وتهع۲ 1 (ه) واه مه ممتاععز ععمتن 66ص بان ۵۲ اممصصه۲ مه ماج صمتاصاجه‌صهع-190 50۵ 1۳6 -17 عباه۲1 
2۲ 


۲۳ نشریه آب و خاک. جلد ۳۴. شماره ۲. خرداد - تیر ۱۳۹۹ 


پس از اعمال شرایط فوق‌الذکر در مدل ۳۸۷۷۸۲ به‌منظور 
بررسی تأثیر سیاست فعلی روند برداشت از آبخوان بر نحوه پیشروی 
شوری مدل برای ۶ سال آینده با گام زمانی سالیانه (مهر ۱۳۹۳ تا 
شهریور ۱۴۰۰) با فرض ثابت ماندن استرس‌های وارد بر سیستم اجرا 
و نتایج آن در مقطع عرضی ۸-۸ استخراج گردید. لازم به ذکر است در 
تمامی اشکال ارائه‌شده محور افقی برابر ۱۵۰۰۳0 و محور عمودی 
7-7 است و دریا در سمت راست واقع گردیده است. همچنین 
خط هم غلظت ۵۰ درصد شوری به‌عنوان مرز مشترک بین آب‌شور و 
شیرین در نظر گرفته شد. 

شکل‌های ۱۳ و ۱۵ به ترتیب روند پیشروی شوری و حرکت خط 
هم‌غلطت ۵۰ درصد شوری برای پارامتر 1195 را در طی زمان نشان 
می‌دهد. نتایج حاصل از پیش‌بینی مدل برای پارامتر 1195 نشان 
می‌دهد که در سال اول پیش‌بینی مقدار پیشروی خط هم‌غلظت 

۰ درصد شوری به میزان ۲۴۰ متر به‌دست‌آمده است. در سال 

سوم پیش‌بینی» میزان پیشروی ۲۶۵ متر شده است که بیانگر افزایش 
۰ درصدی پیشروی شوری نسبت به ابتدای دوره پیش‌بینی می‌باشد 
و در نهایت در انتهای دوره پیش‌بینی میزان پیشروی به مقدار ۳۴۰ 
متر است که نسبت به ابتدای دوره پیش‌بینی حدود ۸۰ درصد افزایش 
نفوذ شوری را شاهد خواهیم بود. 

شکل‌های ۱۶ و ۱۷ نشان‌دهنده نحوه تهاجم شوری و حرکت خط 
هم غلظت ۵۰ درصد کلراید به سمت آبخوان است. همان‌طور که 
مشاهده می‌شود روند تهاجم شوری غلظت کلراید در طی زمان مشابه 
غلظت 15 است به نحوی که در سال سوم و ششم پیش‌بینی خط 
هم‌غلطت شوری نسبت به ابتدای پیش‌بینی به ترتیب به میزان ۲۵ 
متر و ۲۰۰ متر بیشتر در داخل آبخوان نفوذ می‌کند. این میزان نفوذ 
بیشتر ناشی از بهره‌برداری بدون برنامه‌ریزی بوده و سبب هجوم 
جبهه‌های آب شور از طرف دریا به ساحل می‌شود. 


نتیجه گیری 

بحث بهره‌برداری مناسب با کمترین اثرات سوء یکی از مهمترین 
زیرزمینی و شورشدن آبخوان‌ها بخصوص در آبخوان‌های شمالی که 
در معرض تهاجم جبهه‌های آب شور قرار دارند بسیار مهم و حائز 
اهمیت است. در پژوهش حاضر آبخوان دشت ساری-نکا با استفاده از 
مدل عددی ۷101۳101۷ و 11317315 شبیه‌سازی جریان آب 


منابع 


زیرزمینی و انتقال شوری انجام گرفت. پس از ساخت و اجرای مدل 
جریان آب زیرزمینی. مدل کیفی آبخوان جهت شبیه‌سازی انتقال 
شوب اسفاه ارشاا فصو واه ها حرش سل 
5 تهیه و در همان دوره شبیه‌سازی مدل جربان آب 
زیرزمینی و با گام زمانی ۶ ماهه اجرا گردیده و پس از آن پارامترهای 
موّثر بر پراکندگی واسنجی شد. نتایج مدل کمی و کیفی بیانگر دقت 
بالای مدل شبیه‌سازی‌شده است. در ادامه جهت بررسی هجوم جبهه 
آب‌شور به آبخوان دشت ساری-نکا با استفاده از مدل 5۳۸۵۷۷۸۲ 
به‌جای در نظر گرفتن تمامی محدوده مطالعاتی» مناطقی که دارای 
5 با غلظت بیشتر از ۲۰۰۰ میلی گرم بر لیتر بوده‌اند در نظر گرفته 
شد. از مدل ساخته شده در محیط نرم‌اف زار 572۸17۸ برای 
پیش‌بینی وضعیت پیشروی آب شور در ۶ سال آینده با فرض استمرار 
روند کنونی در تغذیه و تخلیه آبخوان. استفاده گردید. تهاجم شوری 
برای غلظت 115 و کلراید رفتاری مشابه داشته است به طوری که در 
سال سوم و ششم پیش‌بینی میزان تهاجم شوری به ترتیب حدود ۱۰ 
درصد و ۸۰ درصد نسبت به سال اول پیش‌بینی افزایش یافت. 
همچنین حرکت مرز مشترک آب شور و شیرین به سمت آبخوان از 
۰ متر در ابتدای شبیه‌سازی به حدود ۴۴۰ متر در انتهای 
شبیه‌سازی رسید. نتایج پیش‌بینی مدل حاکی از آن است که با ادامه 
روند کنونی آبخوان و عدم‌تغییر در شرایط فعلی حاکم و بدون در نظر 
گرفتن افزایش یا کاهش برداشت از آبخوان. امکان بهم خوردن تعادل 
بین آب‌شور و شیرین» پیشروی بیشتر آب‌شور به سمت آبخوان دشت 
وجود دارد. نتایج به دست آمده از تهاجم جبهه‌های آب شور دریابه 
ظرف آیشوای ماخ دهد که این موخنرع ای مطالعات سای 
محققین است. نقش بهره‌برداری و افزایش آن در پیشروی جبهه‌های 
آب شور تأثیر زیادی داشته که نتایج این مطالعه با نتایج استین و 
همکاران در سال ۲۰۱۹ (۲۵) که به بررسی و نقش تهاجم و اینترفاز 
آب شور و شیرین پرداخته بود همخوانی دارد. حفر چاه‌های غیرمجاز و 
اضاقه تردافی از جاههای ابش ترنن غانلی اس کنهشایم آب 
شیرین زیرزمینی در دشت ساری-نکا را مورد تهدید قرار می‌دهد و 
موجب افت تراز آب زیرزمینی و به‌تبع آن تسریع روند پیشروی آب 
شور خواهد شد. با توجه به اینکه با پیشروی آب شور در آبخوان‌های 
ساحلی برگشت آن زمان‌بر و مشکل است لذا لزوم مدیریت و 
تصمیمات صحیح به‌منظور جلوگیری از اثرات مخرب زیست‌محیطی و 
نابودی آبخوان ضرورت می‌یابد. 


6 ۱۷۵۱6 /0۲۵۵۵) 0۶ وععاعصصهتعظ متصصهط و۳۱۲0 ۵۶ ممتامصصت۴۵ .2015 ۷۲ مهتطه‌طه۵۱ ر.ل م0 601 -1 
281-۰ :(19072 ۹616866 رم مضه عماه ۷۷ ۵۶ متام ۱۷۲۵۲۱۱۳۱۱۵۱۷۰ عملونا ۷۷2/۵۲560 56221 -طووهم‌طنامک 1 


(امجتطوهاه طفتاعص طز۳ صوتوعظ م) 


شبیه‌سازی پیشروی آب‌شور در آبخوان‌های ساحلی ۲۸۴ 


طویامطا فاصممتهد‌امومی ۵۶ ممصصعمصظ. فطع متقحته ما فاعم فصو .1979 هط تفص همه ر.ظا جموعلدخض 
97-6۰ :(902 ۰160801027 لمح ممصمزمه له)مممصمم ۱ص ما ومع لممتاتت .عماوره سم تمه رامع 

,نام ۷«حصعاطاد۱۷۷2) موه معلنامه لقلمقاع ویامعصمعمتعامط ه صا )0۲حفصه مصق:مل عصتاهانامطزه. .2010 بک1 أبعلهزصهظ 
۰ 2004 ۵0615 ماصملنه۳ه2 بام‌تاهاتام عصتوتا رصفعتطمز]۷ 

ععج ۷ 921 )عصمهجممعن۱ -افمه۱۳ هم مصتلم0م۷۱ لمممزوممصصزن-عع۲ط1. ,2001 ,2 صفطن) عصمط مه ری مهن عععط 
:(39)1 0۳۵00۵6 فص رعمص۷ ۳:۵ عممصفطه. رصفتوظ من مطقل مد فمتتنامه همهم ۱۷۲۵۱۵۱۱2۵۲۵0 مد ممتفنام8] 
137-3۰ 

مادم طفعتدظ له ماع ما ممتهااصته )عمممصحا مستامی لهج 100 ماج لصا0زن0 :2015 هر مصقو 20 ر.یظ تقلمطان 
0۰ (21)2 همهم متومامعل و ۵۶ تممامز رتاو مقهت تصمصصیم) بصتقام 

جمتفتصا تمه سجعو گم ممتاح‌ابصته آقومتعصمصصتل مععط]. ,2012 ,۲.۳ ۱۷۲۵02 200 ره صهع۵ن1 ریک ما۲ ,۱۷ دام 
262-0۰ :464 ۳۱۷۵۲۵۱۵2۷ ۵۶ امصتتاه0ز بمتهاظ عته)اه(1 بای‌امی) فطا صا پرل‌بتاد عقهع خر :فتعلننامه لقاقدهع 

ممتفاحا مهو ۵ عصتامعتاهع۲۸۷ .2017 بع1 هزیاما۳22 24 رب تلفصطه توطماها2 ره تنه‌مفصطفطه ریم 2220 نطا۲۵ 
آمسته۲ فمد مسانه‌تهه ماطمصتمامنگ مه صیتومم‌صوگ تقممتاجل 6۳ ٩۸۳.‏ عمج آمعتتمصنه عطا عصنون مفتنه عقاه۲ 
(صم1:عظ م1) ,صحیطاه 1 رععمتتاموم۴ 

عمط گم صمتاهاصصه عم همم عنام ه :۳۴۸۷۷۸۵۸۲ مه ع4تاع وتمونا .2002 .نی طابععصقا مه .۷۷ میات 
06-۸7(۰ ۱0۰) ۲10۷ ۲عاه-صامتع متفجمل-۱۵01۱6 ۷۵۲ 2ط10فصمصصتل 

٩2۵1 ۵۲‏ مج م۱ مدب مصنامتع تمممتوع: ۶ ماهنت .2000 ۷۱2۰ مامععوظ هه .و ۷۷۵۵۵ ,0۷۵۵۵96۷ 
772-۰ :(38)5 0۳0۵۱9۱۵۱۷۵۵۵۲ :۳۱۵۲102 تاصنامی مصقصمع۲۱ ما ممتعفتاط1 

عاتصت منود عگننامه لماعدمی ما ممتفااصا ۵20۵و ۵۶ عصتاع0 مه آهعتتمصصنیا ۲ :2006 .و 220 ۷ طمعاطاهه 24 ری تتصصهز۴ 
۰ :(5)6 وعتامصصمطاه] ۷ ۵ فصمت)عععصه ] ممطو ۷۷ .00طمو صفمجهر عناوم مصصتتام۲۷ 

۷۵ 0601081021 ۲۲.9۰ عطا ,۱۷۲۵۲۱۳۲۸۵۷۷2000 .2000 ۱۷۲.۰ ۱۷۲۵120۴۵10 مضه ).۱۷۲ ۳۱۱۱۱ ریکا.۲ هاصظ .۵.۷۷ طویهه ت۲2 
۰ ۲۱۵۱۷۷ 206۲ 0۳۵0-۱۷۱۷ فطا مه فامعمهمن) ۱۷۵0012۲122610 مه 46تدای) بعوت]-۵061 عماج ۱۷۷-صنا۵عن0) تقاتال۷]۱۵ 
٩۰ 06601021021 5۱۲۷۵۲ )92( 4‏ .۲۲ ,۵۵0۲۲ ع11-صو0۵۵ 

-511 طا ممصهاهه ماه 0صنمع ۵ ممتامصتططعاع لمح ممتامتال۳۷ .2009 به طهططافه مه با کنامم‌نج‌طصمدت رب ات2 
بطق ۵۶ عصنمه‌عصنعمظ فص مععصمنه؟ اصمصهعممد]۱ 0عطعهاه ۷۷ چم معصعته‌کجم لقصمنادا 5۳ .صنقاظ 0مصتلصمعصنا ماعل 
(صفتوع۲ ط1) .ت1۵ 

وتعنتامظ )تعوعنا ما ماه )521 ۵۶ ممهم۲ع]۲ع)۱ ۵۶ ممل/2ع 18۷6۵۱ .2017 ۱۷۲۳۲۰ 10صامطتصحظ مضه ,بط صهلهمطوم]۱۷ صهل‌تهکز 
داتعم از صقتوع۴ م1) :673-688 :(3103 5011 20 بعا2 ۷۷ ۵۶ لمصبمز .همم مهو رصعهعه‌مطکا مامو :ناه مفهن)) 
(2090201 

1681014 0601212 طمصصو هو ۵۶ لهتاصهامم وعمتاموع؟ مماهصتامتع 01 ممتاه‌تل۳۷2 .2005 ۱۷۲۰۱۷۲۰ ۸۵۶۵1 یط 011 ر.ط امک 
٩0. ۱۷۳۹] ,-01-۰‏ ]5۵۵0۲ ,1۵001۲۵۲01 فصمت)ه بای لم)مممممز ۲ص م01مصصتا ]۷ 

ص ممتمتاطصا تماه هه وم معصفطل ماقصصئاه ۵۲ ماهع]1ظ .2014 ۰ظ تصوزطاوهعتهاه 20 ,.ن طم0ع۵02مصصطه]۱۷ بیط هام 
(صفتعیع۳ ص) ۲۵۵22۰ ,تمد اه واتعهندتا بععصمتعکمم متام بو صمتصمن؟ 13۳ .عمکننوه آمامدمی نوماه 

فص عصنا ۱۷۲۵0۵ مه میت .فه‌گننامه لماقدمم ووعتعهصا ماه هه صا تماه نامع گم اصمجصم‌ممصه]۷ .2006 0۴۰ عقمصییک 
540-60۰ )صعجصمعه۷]2۵ 

تمانالم۱۷۵ بعتیاو لهعلوما۵ع0) وتا عمط ,۱۷۲۵۲۸۲۲۸۵۷۷2۵0۵ .2003 ۷۷۰ مان 20 ,.ظ ۷۷ تملممجه‌مطو رنلات طابمعمم [ 
5 ۳۱۵۷ 2712016-1(60510۷ ۲۷ مظ) اد م810ع۷ ۳۴۱۸۱۷۷۸۵۸۲2۵0۵0 عظ) ۵۶ صمتاه)صمصتاه ۵061-10 م۱۷۷2 -0صا0۵) 
2003-426(۰ ۲6۰) (۷۲1) ۳۲۵6۵85 ۰۲۲۵60۵016 ۱۲1312۳۷۲5 10۱۵21۵060 عطا 0صج رن ۷) 

۷ ۱۳6 0۶۴ ملظ وعمتناموع ما۷2 مطا ۵۶ )مک رصههلمه2ه۱۷ ۵۶ دمحم ماه ۷۷ لمطمتومک ,برععص ۵۶ رتوتط۱۷ 
(912ع۳ 1) ,2014 وفته‌عصتعص عصتاانفومی 9011 0ج ماج ۷۷ مهتم 1 بر1503 600) 5271-06162 2162 

۵ ۱06۱17ظ فط ما صمتق اما تمحسهعه عصتاام۷]۵0 .2007 ۲ کل 19000ظ مضه یت 9۵1۵۵۵۲ ره یک ۱۱2۲2۷28 
217-۰ :(89)3 )ممصصمهها ماج ۷۲۷ لجعبانه‌تتوخه .همتلهتاوم رصهافممعبم) طاعملظ یحعیه 1۳11826109 

را 0 ۱ ( ۱ 
۰ 0۶ م۲210 ]وع7 اب9۵ تصمنوع۲ لهاعهمع طاً حصعور٩‏ 

(صفتویعظ ص) .عبمکننوه صا عصنام0مظ )ع0حقصعا اصمصتصهاومع مج 10 0عل1مم۸ .1979 ,۸۵ 10002۷1 صه ریک 22۵1 

عاطما ۲ع۷/2 ۵۶ عصمنانلجمه طاً صمتفتتاط ماه امد ۵ وله لممتیمصصتط .2016 ۲۰ تصحبهتطوه هه یه عا۵0 1۵2۵ 
(صفهتوع۲ 12) .صقتطه 1[ بعومرم]ع100ظ مه «عتمص رتصمجصممصظ وم ممممعع]صمن قمع 191 16 0120۷10۷۷7۵۰ 

لماعدهه مطا ما ممتویاتاصا تماج لو عمط ۵۶ ممتاندمم اصمتوصهت مطا گم ممتاماندصزو 1980 ۱۷۷۰ تممملووع5]0 204 ,رت 1۴۵۷۷۵ 
۵۵۵۵۵۰ 216۲ ۷۷ طا و)ممصماط ماتصل اقهمن ۱۷۵0۲۵8 عوعط تعکنبام20 

۲ 006۱ )حرمصماه مفتصن امطمتقصمصصتل-۱۳۷۵ ۸ .1988 ۲۲ تمه اقطهه1ع0طه مصه ر.ظ طفصته دز ۷ ۷ ما۱۷ گلتعط 
.11-۰ :103.1-2 ۳۱۷۵۲۵۱0۵2۷ ۵۶ امصیدام .فم‌گننمه لهاقهمم ط (ناط21۳-۳۲۳) ممتوهمول 

۲ ۵)6/صنا0ع مصلهه عمتم‌صتم ۵۶ ا0ع11م ۰186 .2019 .0 صهلز. 20 یک تمطفهک ریظ «ملقطه ررلا تافنطهع۷ ریق صلهماه 
46-7۰ :156 موم عماج ۷۷ .ومتصصحصل ممهعاص1 ما72 مصتلهف دوع مطا ره مم)مصنلموع 

۵۵۲ 08۴ 16۷61 07۷2۵6۲ص0ع لصح 16۷۵1 962 ۵۶ معصفطله فطا عصتامع 1۳۷6۵۱ ,2013 ۳۰ تصععمطعاهط۸ 220 .۲ ۱۷21201 
.ححتطع۲ بعومع]م100 0هه رععحظ متصمصممححظ ون عممعتهلمم آفممناها 1۳ م۲ .متع‌کننوج 0مصجمعص ما صمتفتاطز 
(۳6۲۹120 

۲ 0۴ هنانز .ممتفااص مهو ۵۶ آمتاجمی تمصصتامن ۲۵۶ [ع0مصظ مصتصصعاظ .1988 ۵ «مصصز۴ 20ظ هه ریک ۷۷11115 
163-۰ :114.2 )مجممههه۷]۱ 0ص عصتصصهاظ وعمعباهوم۴ 


۷۷۱۱۵۷ :۷۵۲۱ ۱۵ .(353 .0 و2 ,[۷۵) عمط )0۲مقصهتا اصحصتحصصقهادمع ۸۵۵1160 .2002 0.1۰ تامصصوظ مضه بت عصعطم -28 
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۱۱۱۱۵۰ 

۶ ۱۵06۱ ]0۲مفصه ‏ ومزمعموتآننمه آمصمتفصمصصیل-ععطا تقلنال‌مصظر 2 :۷۲13121۷8 .1999 .ظ.ظ مصه/۷۷ 200 بی) عصعط2 -29 
مدرد راهطا فاصحصتصصهادم .۵۶ .. فصمتامهع: . لحعتصمطه . مه رهملوهمونه . رعمتامع20۷7 ۵ طمت)ق1 ناه 
۰ ۲ مصصو‌طاهاخ .ع0دناع وتهویه مه ممت)ه)صمصنهم 


نشریه آب و خاک (علوم و صنایع کشاورزی) [م6 0ص ۲۷۵۵ ۵۲ امصنته[ 
جلد ۰۲۴ شماره ۲. خرداد - تیر ۱۳۹۹ ص. ۲۶۹-۲۸۶ | 269-6 0۰ ,2020 ۷۲۵-126 ر2 0۰ ,34 ,۷۵ 


م۸ اهاعدم) طا ممتعصمصا ماو جع ۵ ممتداسصصز 
(وعه صهزم‌عمت 6 ۵۶ معامصا محمطنمو عمط :رتاو ععمت) 


"ص00 هطع۱۷]۵ «مل 1۵ ۲۲۱ -* ز0تصها! ۱۲۱ - نزمه :۸۲ 
15-04-19 :16061۷60 
09-02-2020 :666060 


فص عمتهلتعجمی لامج لهاصمصصهمتجه. مه ممجمععه فقط صمتفتااصا مج ۵2و ۵ متاووز م1 :ممتام00 1۳۸0۲۵ 
۶ 1201 0ص ممتامازماوه همع 60وهع۲مه۱ ,فم‌گننامه مادم صم ممتامهتان ماه 0صتامتع ما فصع معصلفهع1001 
06 10 160 عقط ممتفتصا مظ1 تمکننمه لماودمه ماصاً ممتعام ماه ۵2و 60فننهه 276 عصمجم‌معفصهمظ فععتناوه لقافدهه 
فقط عمتورصتم فنطا .عمما0و۲۵ م۳2 0۴ حمتاعااماه فطع ما فصصفامامم تصقوه عصتفنی . رفته‌گتنامه 0 صملا۹21112 
۵71-۵: 6 صا متلمان اهوم ۵ ممتاه)متهاع0 مه ممتفتاطاطا مله ۵2و 0۶ کلف مط 0عووعمص «القتمتاومع 
۷۰ ۱/۵۲6۲ 16ات۵۱ مه رد10 رمتناتام‌آوه 1۵۲ ماه صامع ۵۶ معموه عطع فاتحصتاً ممتقاصا تمه سدع .ما20 
ص 106200 و1 معت2 تاد مج[ .0ع0ع۵۱و مه تمئننامه م1م-۹271 ,رتاو اجعومزم مط ۳۵ :۱۷۲۵0۳۵09 0ص2 واهزه)۱۷2 
0 ۱560 15 2000 ط۷6۵10 ۷۲۵۱۳۲۲۸۵۷۷ م1 .صویا ۵۶ ابهم مدمه فطا ط بهعو صفتم‌فمی عطا ۵۶ فععمطه صتعطانامد فطا 
,2060016 نامز متقله 10 عنام م)اه- 21 ما موه اما عماج 0صنامنع اصمتفمه) 4صه مه)ه-۹0620۷ 2 ع)هاناحطله 
مه علصتومیامی مه رصملوهمونه همهم مملامع07ه ,۷۲۲31۷۲۹ ما .عون و27 5۳۴۸۷۷۸۵۸۲ 20 1۷۲1312/15 
]۳۸۷۷۵ رفتاط1 .10۷عجع عط) )عع1ه مطا عع0نقصمن )مه فعمل ۱۷۲۲312۳۷۲۵ 6۵06 ل2مزععصصتنه مظ1. ,0عفا وت2 220129 
6 92 1107 نامع ۵۶ عصتاممو 10۲ هل ۵۵6982 م1 ,1810211260 فه۳ معقامه ۳۱۵۷ 0ادصونا ماطهتته ۲۷ 
0 ,م2 امدتو۵وها رجمتاهمتمامه هل متصصهطول ملظ مقاقل لهمنع0۲010ع)۵ مقاقل ۲۵80۲۲6۵5 ماع ملصا 011260ععم2ه 
مصمتانل‌همن فانصا ررماقاه )صمتمصهه مصه ماما 20ماو) ۵۱۵۵ ما رامل‌ممهط 100 عطق 10ت0 10 مصمتاممصمکصا لمعنعماهعو 
40 (۲2۵عج۵ع) مصمتاتومن هصتهها ۵۶ ۵۵ مه (ماهای-20عاه معط 10۲ 2010 تعحاصماوعه ما 167۵ مماها2۲0) 
همع حصم فهمتنمم فوعمه ت«لطممصه 48) ممتاهمتامعون ل۲2ممجصعا ر(معقآمهم طظ) مممامدم ۶10۳ ر(ته‌صنمه 
6 0160۵۵۲6 10 ,20061 عط 60 28912060 ۷۵۲۵ ومماد عصصتا زاط‌صمحه مضه (صمت)تلجمع )صمتعصمع صا 2014 ]وناعنه م 2010 
0)) )واه عتهتا50 200۷6 3694 ماصاً رللمه۵۲120ظ 1901601060 و۷۷۵ 2762 «جاه فط رقتیع اع0مصظ )0۲مقصهتا 2 110۷ 
۲ 6 مه ععهتاه تعلننامه ما وم فععصقطه 20۷۵)هتان مط) مققااصته م6 ۱۵۵0 ق۵ 1ع0مهصه ۱۷۲1312/5 مظ]. 0(۰ظ 
۷ 0۲60۵160 ۷۷2۵5 ]020500۲ نامز ۵۶ ام0مه آهتتانمعوهع مط ۲۳۵۲۵10۲8۵ .۲ع1لنامج معط 0۶ طاوع0 مط لته مه 2۵06 
5 مه م06 168۷۵ ماه0صناهع 2010 ععحامجماوعق .1ع8مصه سمل تمتتامهعومع فطع صاً وععصقطه «تدووعهعه مطا عصتلم 
2 024012 16 .1۵506001۷61۲ راع0مظ )0۲صقصه) لمح سمل عطا طا فصمتاتجمی اهتانصا معط عم صمله) مه 0212 01 0هه 
تصعاعومی هه ۱۵۵/1 35000 ۵۶ صمتاح‌اصععوهن 115 صهاوصمی مه ها ۲۵۵۲۵5۵۵6۵0 ور عمط مطا ما عت2 تاد عمط مصتنعلتهه 
,0۵0001۷ 1۲2116 بو رقلاه۳ عمتاهبجعهوهطاه فطع ۵۶ 45فعظ میج عمط 0عتماهه 7۵ راممصظ فتطا م1 .260/1 01282.2 
قانهعامصط 0 اصملمتللهم راتوم۵ رعوعمطهع۲ ععفننامه هه ,هون نامه رولدم۵۲ 611600۷6 راصمام006111 ٩0۲2826‏ 
۰ ۳228۷۵۲۵۵ ۷۵۲۵۵ راز وونل موم عصقا لماه0 ۳۱۵۲2 ,۱۱0 1وممعن اممتهبمتعمم رجمتقتن مه 
صعتعصهه مه ماه طامطا عم ام مدیم وج عمج ما عط ۵۶ ومتله‌طاناهه م1 +صمزوع ور 20 وقلنعع1 
5 0056۲۷۵0108 14 طم هل ۵۶ واهبه۱ میج صمو اطمم امجمتومج من مصه مت مطا عصتوتا فممتاتل‌جمع 
۱ 
بطمتاه 2110 عطتت1۳ .2010 عطجصمامم9 ما فلع۱ تمه نامع که ومباد۷ 0ماقاناملهه همه وام۷ع1 ماه نامع 0۵96۲۷۵0 
۵ ,۱۷۲۱۵۱۳۲۲ 0۷ همه فماهب ع7ع۱ مماه موم انا ودره م۵ مفمنلهب ااامنقجمع مننیه بط 
مط بطمتاه‌طانله م26او- ۱620و ماه رماع پممتهعع ۵۶ 1۵۷۵ ما0هاوهععه مه نصا ومبله۷ ۲۷۵0عوما۵ فطا 2266060 
۵ 02 2014 )وناونض مصه 2010 عطمصماوع م060۵ ۵۵۲100 تقع۷ عنام معط ۶۵۲ 0ماهاتاحطژه و2 [06مقظ )صمتفصه 
مصتااناوع: مطا رعمته‌تطانلت 061مصظر ما ۵۶ لجع فطا ۸ عمماه مصت «اطتصمصه له قمتوج ووعته 48 مفصد 01۷1060 
4 1۲01 ۷۷۵۲۵ ما۷2 ۷۱۵۱0 تمد عصتاانادع۲ مطا ملظ ریجمل/۳ 6 .21 مه 5.3 صمظ 0ممصد تامهم منیه بط 
72) 2014 ۸۵۲05۲ 20 2010 هام56 مهع0۵07 ۵۲00 مطا 10۶ ماما نامه و77۵ [ع0مصظ 10 صمتا۷21102 ۰1۳6 .5.9 ٩2۵‏ 0 
0 224 20001 )ممتوصهتا رام مهای- 0۵207 عط صا کصمتم‌تللهم ممتاهاهمن مطا ۵۶ ومبالغ۷ م1۳ ۰(رو0مته ووعته 
اصعصصعع‌تعه 2000 2 ماجعافتالاً فلانافع؟ م1 .راوبتامع(٩ع۲‏ ,0.97 208 ,0.98 ,0.99 8مصنه)9اه ۷۷۵۲۵ 061مهظ ۶۱0۷ مط صد 2۵061 
بل 220 و11 ممئوت نله فهب (ع0مص )0۲مقصهتا م1 ,واعبه۱ عماج صنمم 0عام1نهلهه هه 7۵0عوها۵ عط م۵0۳۵ 
عطا هصتامءه فعامصصه‌تدم مصتامد ره را (ول‌متنعم ووعتاه 8) 2014 ]5باع۸ 6۵ 2010 ععطاجصمامعه حصم هم ممتام‌تامهع‌صمم 
0 2014 )ناخ حطمظ مق ممتاحتامهعصمع بان 0هه دبا بلعمصظ مطا ۵۶ «ممتنهعج مط) تدم 10 .98ععمعم صمنفتهمول 


01 0۶ م۲۵ روهنهمنااگ متلنیهل ۲ مه مصتمم‌صتمص عماج ۷۷ ۵۶ ۲۳۵۲۵990۲ )صهاووو۸ 20 0۳200۵160 ,۱۷۲,5۵ -2 20 1 
صد ,امصوظ ررعمامصطمع 1 ۵۶ بلومب دنا تمه بتطوواظ امصقظ رع‌صتهعصتم۴۸ 

(01).26.11 6۵ 41 تحصقط :اتقصصظ میمصت طممدمتزمن -۶) 

مه بصحتطه آ رعتمظ۴ ۵۶ رتفتطتا۱ رباعم دممعومط ماه ۷۲۷ رطامتهعومک وممعنامومع عماج ۷۷ 0۶ اصمحص م126۵ -3 
9 :]10 
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عصصم)م5۵ ما صمتاه‌تاصهمهمن 01 هه 11۳5 مطا مصتیهتوممی برظ .عموممتام متام ت۷۵ ۶۵۲ ۵0و ۷/۵۲۵ 2015 همع 
2 طمبامط) 60ع۲عنط20 ق72 ممتاه‌طتلم اع0مصظ )۵۲جقصه1. بصار ۷/۵8 [06محظ عصعتعصه) معط رعمتاتصمی لفتانص عط مه 2010 
فصح 0.83 0مصنه)ط0اه مه 1125 ۶۵۶ 1ع0مظ 0۲۲حقصعتا فط) ما اصمت‌گهمه ممتاقامترمن فط ۵۶ ففتاه۷ 14 .2001-61701- 1121 
معط ما اصمت‌تهمه ممتاقاعمع فطا ۵۶ فمله۷ م1 .1۵5۵600۷۵ رباعم ممتاهلتله۷ مضه [ع0مصظ عصمتفصهع عطا ظا 0.87 
۷۰ ,۱۱0061 20108 تلج۷ صرح 061مه )صعتقفصهع فطع ما 0.86 24 0.82 0مصله)0ا۵ ۱۷۵۲۵ نآ ۶۵۲ 061مقظ )28201 

2 مصتجصه۲ عممتاتصمع ملع1۲02۵010 7 مط له ممتصصاووج مه ۵061 )تمموقصه ۵۴ ممتاهل‌نله۷ عط ماخ خصمتعتمومن 
لماعدهه عط مفصا صمقتتاصا تماه 62و تعطاعبظ اقط) وماهع0ه1 عم ۳۸۱۷۷۸۵۸۲ مصتوه صنا ممتاهانه کقهبد6 2۲60100۷۵ 
1 0۷و عصا صد کتامعن صهع فععلنا20 


جماص مه 9۵20۷ رعلتنامه م)له اه رصمتاه‌اصهعصمع ]921 رتمکنامه آماعهمن 1۵۲۱۷۵۲۵۹۰ 


